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We study the effect of dilute pinning on the jamming transition. Pinning reduces the average 
contact number needed to jam unpinned particles and shifts the jamming threshold to lower densities, 
leading to a pinning susceptibility, Xp- Gur main results are that this susceptibility obeys scaling 
form and diverges in the thermodynamic limit as Xp oc |0 — where cj)^ is the jamming 

threshold in the absence of pins. Finite-size scaling arguments yield these values with associated 
statistical (systematic) errors Xp — 1-018 zb 0.026(0.291) in d = 2 and Xp — 1-534 zb 0.120(0.822) in 
d = 3. Logarithmic corrections raise the exponent in d = 2 to close to the d = 3 value, although the 
systematic errors are very large. 


In jammed packings of ideal spheres, particles are 
locked into position by their repulsive interactions with 
their neighbors, which in turn are locked into position 
by their neighbors, and so on, so that the entire system 
is mechanically stable. Pinning is an alternate way of 
locking a particle into place, so the interplay of pinning 
and jamming can potentially lead to interesting new be¬ 
havior. Pinning is known to have a rich interplay with 
glassiness; pinning raises the glass transition [T] and can 
be used to probe its nature 013 ] and associated length 
scales 013. In jammed systems, pinning lowers the jam¬ 
ming density 0|6] and allows access to length scales 0. 
Here we show that the addition of quenched disorder in 
the form of random pinning has a singular effect on jam¬ 
ming. In the dilute pinning limit, jamming is highly sus¬ 
ceptible to pinning, with a “pinning susceptibility” that 
diverges at the jamming transition as a power law in the 
thermodynamic limit. 

In spin systems such as the Ising model, the magnetic 
susceptibility can be calculated by considering the re¬ 
sponse to “ghost spins” or especially-designated spins [8] 
in the limit that their density vanishes; similarly, in per¬ 
colation or correlated percolation one can calculate sus¬ 
ceptibilities to “ghost” sites or bonds Bm that are 
vanishingly probable. Here, we consider the response to 
“ghost pins.” Systems of N particles, of which a frac¬ 
tion Uf are pinned, are prepared by quenching infinitely 
rapidly from infinite temperature, T = oo, to T = 0 at 
a volume fraction (f. We calculate the fraction of such 
systems that are jammed, or equivalently, the proba¬ 
bility that a state prepared in such a way is jammed, 
Pj{(j),N,nf). We then define the pinning susceptibility 
in the limit of vanishing pinning: 


Xp = 


dpj{^,N,nf) 
nf^O duf 


( 1 ) 


We find that Xp(0, A") and the probability of being 
jammed, pj{(j),N,nf) obey scaling form and that Xp 
verges in the infinite size limit. 

There have been two distinct approaches to study¬ 
ing scaling near the jamming transition: in terms of a 


configuration-dependent or infinite-system critical point. 
Each finite jammed configuration of particles, A, has its 
own critical volume fraction (j )^, which converges to a sin¬ 
gle value (j)^ only for infinite system sizes [nHl3]. Eor 
many purposes, scaling behavior near jamming is best 
done by measuring the deviation from the configuration- 
dependent critical density (as suggested by Refs. mun]). 
Here, since we are studying the convergence of the distri¬ 
bution of the configuration-dependent critical densities 
to the infinite-system critical density, we naturally make 
use of the other approach, scaling in terms of deviation 
from (j)^. The existence of two distinct scaling pictures 
is seen in many other systems with sharp, global transi¬ 
tions in behavior, as originally discovered in the depin- 
ning of charge-density waves [isHni- Such systems may 
not obey the inequality between the correlation length 
and dimension z/ > 2/d derived for equilibrium systems, 
unless analyzed using deviations from the infinite-system 
critical point [nmn]. 

To study the pinning susceptibility, we conducted 
numerical simulations on packings of N repulsive soft 
spheres in d dimensions mm at fixed area (two di¬ 
mensions, d = 2) or volume (d = 3) in a square (d = 2) 
or cubic (d = 3) box with periodic boundary conditions. 
We considered 50:50 mixtures of disks (d = 2) or spheres 
(d = 3) with a diameter ratio of 1.0:1.4. Particles i and 
j with radii Ri and Rj interact with pairwise repulsions 
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with a = 2 (harmonic) or | (Hertzian) . 

The initial position of each particle was generated ran¬ 
domly, and the positions of Nf = UfN particles (chosen 
at random) were fixed. Gonfigurations in which fixed 
particles overlap were excluded. We then minimized the 
energy of the system to obtain jammed packings. The 
upper inset of Eig. 2 contains a sample configuration for 
N = 256, A/ = 2. 

We calculated the jamming probability Pj(0, N,nf) at 
T = 0 for systems of size A = 600,1000, 2000 and 4000 
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in d = 2, and N = 800,1600,2400 and 3200 in d = 3. 
For the small fractions of fixed particles studied here, 
the criterion for judging a system to be jammed is the 
same as in previous studies [20]: jammed systems have 
positive bulk moduli, energies and pressures. 

As for systems with no pinned particles, we find that 
systems with dilute random pinning are isostatic at the 
onset of jamming. In agreement with earlier studies 
on jammed hard sphere systems with dilute pinning |6], 
we find that dilute pinning can result in a generalized 
isostatic condition. One must distinguish between two 
types of contacts: the number Nmm between two parti¬ 
cles which are mobile during equilibration, and the num¬ 
ber, Nmf, between one mobile and one fixed particle. 
Each of the Nf fixed particles requires no contacts to 
be stable, while each of the Nm = N — Nf mobile par¬ 
ticles require, on average, a minimum of Ziso contacts. 
When Nf = 0, Nm = N and Zi^o = 2d - ^ + 
where the second term arises from d zero modes asso¬ 
ciated with global translations allowed by translational 
invariance, and the third term is needed to ensure a 
nonzero bulk modulus [20j. Our equilibration protocol 
breaks translational invariance when Nf > 1; thus in this 
case Ziso = 2d ^. Since the average number of con¬ 
tacts for a mobile particle is Zm = {2Nmm + 
the number of excess contacts that constrain mobile par¬ 
ticles is N excess = Nm{Zm “ Or 

^excess — ^ram ^mf dNjyi dQ 1 (3) 

where q = 1 ioi Nf = 0, and q = 0 for Nf > 0. Fig. la 
shows that this relationship is upheld: isostaticity means 
that the excess number of contacts approaches zero as 
p ^ 0. Additionally, Fig. la shows scaling collapse onto 
universal curves as function of rescaled pressure, p^^‘^N. 
This is exactly the same as what is observed for systems 
without pins m- Note that Fig. la is analogous to Fig. 
2c of Ref. [20| in the absence of pins, which shows not 
only the region of slope 1, but crossover to a slope of 2 
at very low pressures, arising from a Taylor expansion 
of {Z — Ziso) in P foi* finite systems. In Fig. la there 
is perhaps the hint of a crossover to a higher slope at 
the data are quite noisy at such low 

pressures. 

One might expect that since pinned particles support 
mobile ones, the number of mobile contacts will decrease 
with increasing pin density. Indeed, for all pressures 
studied, increasing Uf decreases the average number of 
mobile contacts, (The only exception is a slight, 

but nevertheless reproducible, uptick between Nf = 0,1, 
related to the loss of translational invariance.) Fig. lb 
shows Zr^ broken down into Zrnm arising from mobile- 
mobile contacts, Zmm = 2Nmm/Nm^ and Zmf from 
mobile-fixed contacts. At a pressure low enough to ap¬ 
proximate the jamming threshold (circles), the average 
contact number Z^(n/) (filled circles) is well-fit by the 
dashed line shown. Interestingly, raising the pressure by 
a couple of orders of magnitude does not result in sig¬ 
nificant changes in Zmf{Tif) (red symbols). In Fig. lb. 






FIG. 1. a: Excess contacts as a function of rescaled pres¬ 
sure. d = 2 and d = 3 systems are as marked. N = 
600,1000, 2000,4000 colors are red, green, blue, black. N = 
800,1600, 2400, 3200 colors are red, green, blue, black. Sym¬ 
bols for Nf = 1,2, 3,4 are square, circle, triangle, and dia¬ 
mond. b: Contributions to contact number for mobile par¬ 
ticles. Zmm — 4.0 are shown as blue, Zmf as red. Symbols 
for logio p = —6, —5, —4 are circles, triangles, squares. Filled 
black circles fit by dashed line are Zm — 4.0 for logio p = —6 . 

Zmm — 4.0 is contrasted with Zmf to show that mobile- 
mobile contacts disappear more rapidly than mobile-fixed 
contacts replace them. Thus, jamming in the presence of 
pinned particles is an unexpectedly “frugal” process, in 
terms of its use of mobile particles to produce global me¬ 
chanical stability. 

Increasing Uf raises the probability of jamming at any 
given value of 0, in accord with previous work on jam¬ 
ming in the presence of fixed particles lai- Increasing 
N steepens the jamming probability, as in the absence 
of pinning m- These features are illustrated in Fig. 2, 
which shows the jamming probability pj{(j)^N^rif) aver¬ 
aged over 10,000-30,000 d = 2 systems of size N = 600 
and 2000, with harmonic repulsions and Nf = 1,2,3 and 
4 fixed particles. The dashed lines in Fig. 2 are fits to a 
two-parameter logistic sigmoidal form: 

Pi ( 4 ) 
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where a{N,nf) is the “width” of pj^ in that it spans 
probabilities from ^ to |; while b{N,nf) is the value of 
volume fraction at which pj = For all N and Uf 
studied, logistic sigmoid fits to 13 — 21 independent (j) 
values result in values less than 0.10. A slightly more 
flexible three-parameter fit, to “Richard’s curve”, does 
not yield significantly better measures of goodness-of-fit. 


limit n/ —> 0, or ^ ^ 0. (We note that in fitting these 
quantities to our numerical data on pj vs. 0, the limit 
Uf —> 0 is taken as the limit Nf 1 and not Nf 0 
since one pin destroys translational invariance. The dis¬ 
tinction between Nf = 0 and Nf = 1 is irrelevant in the 
limit N ^ oc.) We seek the behavior of d{y) and b{y) 
near ^ = 0: 



0 


FIG. 2. Probability pj of jamming as a function of packing 
fraction, 0, in d = 2 systems for N = 600 (red) and N = 
2000 (blue). Nf = 1,2, 3,4 are represented by square, circle, 
triangle, and diamond symbols, respectively. Dashed curves 
through data are fits to a logistic sigmoid function. Upper 
inset: An equilibrated configuration in d = 2 with Nf — 2. 
Lower inset: Value of 0 such that pj — 1/2 ior N — 600 (red) 
and 2000 (blue) versus A/, the number of pinned particles. 

We now propose a scaling ansatz for pj . Since the frac¬ 
tion of pinned sites, Uf, is an independent control param¬ 
eter with which to approach 0^, a two-variable scaling 
function can be constructed for the jamming probability. 
There is significant evidence that the upper critical di¬ 
mension of the jamming transition is dc = 2 [201423] . For 
d > 2, we therefore expect that finite scaling depends not 
on linear system size, L, but on particle number N [24] . 
We therefore propose 


Pj =F{N(t)N^,nfN^^^) ( 5 ) 

where A0 is the distance from the jamming transition 
for the unpinned, infinite system: A0 = ^ — 0^- 

We can rewrite Eq. in terms of the scaling variable 
X = N(j) N^ as 


Pj = 


1 


(1 -h 


( 6 ) 


with d = a and b = (j)^ — ab . From the logistic sig¬ 
moid fits we can obtain functions d(^), b{y), critical ex¬ 
ponents 7 p, and the jamming threshold 0^, where 
y = UfN^p'^. Because the pinning susceptibility is de¬ 
fined in the dilute pinning limit, we are interested in the 


d = aoFaiy] b = bo F biy (7) 

with higher-order terms in y neglected. Table I shows 
the parameters in pj from nonlinear least square fitting 
in d = 2 (four system sizes, four pinning densities) and 
d = 3 (four system sizes, two pinning densities). 



d = 2 value 

error® 

error^ 

d = 3 value 

error® 

error^ 

V 

0.491 

0.004 

0.045 

0.439 

0.007 

0.051 

Iv 

1.018 

0.026 

0.291 

1.534 

0.120 

0.822 


0.8419 

<0.0001 

0.0001 

0.6472 

0.0001 

0.0004 

ao 

38.555 

1.088 

12.329 

46.177 

2.558 

17.521 

bo 

1.648 

0.017 

0.193 

3.370 

0.065 

0.442 

bi 

9.646 

0.879 

9.958 

2.970 

1.145 

7.844 


TABLE 1. Best fit parameters for Eqs. [^ “Traditional 
statistical errors. ^Bounds roughly incorporating systematic 
errors. Estimated from deviations between model and data 
m, systematic variance is approximated as twice the best-fit 
divided by the number of parameters. 


The parameter ai is sufficiently close to zero that it 
is not listed in Table 1. The widths of the sigmoids do 
not vary significantly with Nf] the principal result of 
increasing Nf is move the sigmoid to the left, leading to 
jamming at lower values of 0 (as in lower inset of Eig. 2). 

The scaling ansatz for the jamming probability, Eq. 
suggests a one-variable scaling ansatz for the pinning sus¬ 
ceptibility. 


Xp(A</.,7V) = |A0|-> QdiAcpN-). 


(8) 


Note that we have an explicit form for the function gd^ 
where d > 2 is the dimensionality, via Eqs. l,[^and[^ 


gd{x) = 


-|- g&o+ao x^2 


(9) 


We also calculate Xp each value of d and N using a 

AT 


Xp 


finite-difference version of Eq. 

for limitingly small values oiN^ Nf. Using = 
2, Nf = 1 in the finite difference yields smooth curves 
for A = 600,1000, 2000 and 4000 in d = 2 in Eig. 3a, 
and for A = 800,1600, 2400 and 3200 in d = 3 in Eig. 
3b. We have additionally calculated Xp using pairs of Nf 
values other than {2,1} (not shown). We find that the 
finite-difference approximation to Xp develops a progres¬ 
sively higher peak as Nf ^ Nf. Uncertainty in XpW 
arises from error in the parameter a (error in b con¬ 
tributes much less) which is fit independently for Pj (0, 2) 
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FIG. 3. a: Susceptibility calculated as finite difference for 
d = 2. N = 600,1000, 2000, 2000 with errors as areas around 
curves in red, green, blue, grey. For N = 600, solid line har¬ 
monic repulsion and dashed line Hertzian repulsion; all other 
values of N show harmonic repulsion, b: Susceptibility calcu¬ 
lated as finite difference for d = 3. N = 800,1600, 2400,3200 
in red, green, blue, grey, c: Universal scaling function (curves, 
Eq. and finite-difference approximation (points) for non¬ 
singular part of pinning susceptibility, gd{x), for d = 2, 3. 


and Pj(0,1) before the difference is calculated. This un¬ 
certainty is shown as an “envelope” about each curve in 
Figs. 3a,b. Note that Fig. 3a contains data for both 
harmonic (solid line) and Hertzian (dashed) potentials 
for N = 600. The disagreement between the two curves 
is significantly smaller than the error for either curve, 
supporting the expectation that the pinning susceptibil¬ 
ity near criticality is independent of the details of the 


repulsive potential. 

Eqs. 1^ and [9|arise from differentiating the scaling 
form for pj in EqTl^with respect to its second argument. 
In Fig. 3c, we show the universal functions Qd as curves 
for d = 2,3. Finite-difference results for different N are 
shown as points for selected values of There is 

excellent agreement between the points and the curves, 
indicating that the data at each N are in good agreement 
with the fitted parameter values in Table I obtained by 
fitting to data at all N and Nf. The universal function 
Qd peaks at X = —0.043, —0.078 for d = 2,3, respec¬ 
tively. The scaling form of Eq. implies that in the 
thermodynamic limit, we obtain a power-law divergence 
of Xp ~ |A</)|“>. 

Note from Table I that the values of (j)^ are in ex¬ 
cellent agreement with previous work on bidisperse soft 
spheres mun]. The finite-size exponent v ioi d = 2 
and d = 3 is consistent within uncertainty, as expected 
for systems at or above the upper critical dimension. It 
is also consistent with the central-limit-theorem value of 
V = 1/2, identified earlier in the absence of pinning in 
Ref. [TT], and with v = 0.465 ± 0.01 obtained for d = 2 
systems by Vagberg, et al. m, who included power-law 
corrections to scaling in their analysis. 

The pinning susceptibility exponent 7 ^ in d = 2 and 
d = 3 is significantly different when we consider only sta¬ 
tistical errors - in contrast to the dimension-independent 
values expected above the upper critical dimension of 
two. However, these do not include systematic errors 
due, say, to choice of theoretical analyses. For example, 
one way of including logarithmic corrections in d = 2 , 
the expected upper critical dimension, leads to a consid¬ 
erably higher value of 7 ^ = 1.50±0.95 (statistical errors). 
This agrees well with 7 ^ ^ 1.53 in d = 3, but one must 
recognize that the estimated range of systematic errors is 
enormous. As a proxy for exploring different theoretical 
models. Ref. [25] proposed a method which explores the 
range of fits that is comparable in residual to the best 
fit. Following their prescription, we find much larger sys¬ 
tematic uncertainties in our estimates of 7 ^ (Table I). 
Therefore, our numerical results cannot resolve whether 
jp is the same in d = 2, 3. 

Indeed, one can argue that 7 ^ may depend on d as well 
as a d-independent exponent, u. Conceptually, the nar¬ 
rowing of the jamming transition with increasing system 
size m and the shift in the average transition 016], 
imply a derivative of the jamming probability which de¬ 
pends singularly on the density of frozen particles. For 
attractive pins in d = 2 [5], it was noted that average 
distance between pins. If could be equated with a cor¬ 
relation length ^ oc at the jamming threshold. 

Since If oc this argument suggests Ac/) oc 

as a scaling relation. Our Eq. would thus be writ¬ 
ten Pj = F{A(j) N'^ implying 7 ^ = du. With 

= 1/2 in d = 2, 3, this d-dependent prediction for 7 ^ is 
consistent with numerical results of Table 1. 

In summary, we have found that jamming is infinitely 
susceptible to pinning at the jamming transition in the 
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thermodynamic limit. We have identified a new expo¬ 
nent associated with power-law divergence of this pinning 
susceptibility. The divergent response to pinning, even 
in the limit of infinitely dilute pinning, suggests that it 
should be fruitful to study the interplay of jamming and 
pinning at higher pinning fractions. 

We thank A. A. Middleton for the scaling argument 
leading to 7 ^ = dz/, and thank R. Kenna for instructive 
discussions. Acknowledgement is made to the donors 
of the American Chemical Society Petroleum Research 
Fund for support of this research, and to the Division 


of Natural Sciences and Provost’s Office of Swarthmore 
College (ALG, SN, EP). This research was also supported 
by the US Department of Energy, Office of Basic Energy 
Sciences, Division of Materials Sciences and Engineer¬ 
ing under Award DE-EG02-05ER46199 (AJL, GPG), and 
by the National Science Eoundation under award DMR 
1312160 (JPS). This work was partially supported by a 
Simons Investigator award from the Simons Eoundation 
to AJL and by a University of Pennsylvania SAS Disser¬ 
tation Award to GPG. ALG was the recipient of a Eugene 
M. Lang Eaculty Eellowship from Swarthmore Gollege in 
support of sabbatical leave. 


[1] K. Kim, Europhys. Lett. 61, 790 (2003). 

[2] C. Cammarota and G. Biroli, PNAS 109 , 8850 (2012). 

[3] L. Berthier and W. Kob, Phys. Rev. E 85, 011102 (2012). 
[41 S. Karmakar, E. Lerner and 1. Procaccia, Physica A 391 , 

1001 ( 2012 ). 

[5] C.J. Olson Reichhardt, E. Groopman, Z. Nussinov and 
C. Reichhardt, Phys. Rev. E 86, 061301 (2012). 

[6] C. Brito, G. Parisi and F. Zamponi, Soft Matter 9 , 8540 
(2013). 

[7] M. Mailman and B. Chakraborty, J. St at. Mech. 2012, 
1201.0812 (2012). 

[8] R.B. Griffiths, J. Math. Phys 8, 484 (1967). 

[9] P.J. Reynolds, H.E. Stanley and W. Klein, J. Phys. A: 
Math. Gen. 10 , L203 (1977). 

[10] J.M. Schwarz, A.J. Liu and L.Q. Chayes, Europhys. Lett. 
73 , 560 (2006). 

[11] C. S. O’Hern, L. E. Silbert, A. J. Liu and S. R. Nagel, 
Phys. Rev. E 68, 011306 (2003). 

[12] P. Chaudhuri, L. Berthier and S. Sastry, Phys. Rev. Lett. 
104 , 165701 (2010). 

[13] D. Vagberg, D. Valdez-Balderas, M. A. Moore, P. Olsson, 
S. Teitel, Phys. Rev. E 83, 030303 (2011). 

[14] C. S. O’Hern, S. A. Langer, A. J. Liu and S. R. Nagel, 
Phys. Rev. Lett. 88, 075507 (2002). 


[15] C.R. Myers and J.P. Sethna, Phys. Rev. B 47 , 11171 
(1993). 

[16] C.R. Myers and J.P. Sethna, Phys. Rev. B 47 , 11194 
(1993). 

[17] A. A. Middleton and D. S. Fisher, Phys. Rev. B 47 , 3530 
(1993). 

[18] F. Pazmandi, R. T. Scalettar, and G. T. Zimanyi, Phys. 
Rev. Lett. 79 , 5130 (1997). 

[19] J.T. Chayes, L. Chayes, D.S. Fisher and T. Spencer, 
Phys. Rev. Lett. 57 , 2999 (1986). 

[20] C. P. Goodrich, A. J. Liu and S. R. Nagel, Phys. Rev. 
Lett. 109 , 095704 (2012). 

[21] P. Charbonneau, E. 1. Corwin, G. Parisi and F. Zamponi, 
Phys. Rev. Lett. 109 , 205501 (2012). 

[22] P. Charbonneau, J. Kurchan, G. Parisi, P. Urbani and 
F. Zamponi, J. Stat. Mech, P10009 (2014). 

[23] P. Charbonneau, E. 1. Corwin, G. Parisi and F. Zamponi, 
Phys. Rev. Lett. 114 , 125504 (2015). 

[24] K. Binder, M. Nauenberg, V. Privman, and A. P. Young, 
Phys. Rev. B 31 , 1498 (1985). 

[25] J.J. Mortensen, K. Kaasbjerg, S.L. Frederiksen, J.K. 
Norskov, J.P. Sethna and K. W. Jacobsen, Phys. Rev. 
Lett. 95 , 216401 (2005). 



